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Abstract. The long— term, highly periodic and correlated variations in both the pulse shape and the rate of 
slow— down of two isolated pulsars (PSRs) PSR B1828— 11 and PSR B1642— 03 were discovered recently. This 
phenomenon may provide evidence for "free precession" as suggested in the literature. Some authors presented 
various kinds of models to explain this phenomenon within the framework of free precession. Here we present 
an accretion disk model for this precession phenomenon instead. Under reasonable parameters, the observed 
phenomenon can be explained by an isolated pulsar with a fossil disk. This may link radio pulsars and anomalous 
X— ray pulsars (AXPs) and present an indirect evidence for the existence of the fossil disk in nature. 
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1. Introduction 

Two isolated pulsars (PSR B1828-11 & B1642-03) 
present the phenomenon of long— term, highly periodic 
and correlated variations in both the pulse shape and the 
slow— down rate (for PSR B1828— 11, its variations are 
best described as harmonically related sinusoids, with pe- 
riods of approximately 1,000, 500 and 250 days, see Stairs 
et al. 2000; for PSR B1642-03, its timing residuals exhibit 
cyclical changes with amplitude varying from 15 to 80 ms 
and spacing of maxima varying from 3 to 7 yr, see Cordes 
1993 and Shabanova et al. 2001). This phenomenon pro- 
vides the most compelling evidence for precession (Link 
& Epstein 2001). Some authors (Jones & Andersson 2001, 
Link & Epstein 2001, Rezania 2002) proposed different 
models to explain this phenomenon within the framework 
of free precession. However, further investigations are ex- 
pected on whether the force of vortices pinning to the stel- 
lar crust is strong enough to damp out the free precession 
of neutron stars (see the conclusion section for details). 

The pulsating X— ray binary Her X— 1 also shows an 
unusual kind of periodic timing behavior in the form of an 
on— off cycle about 35 days. Sarazin et al. (1980) showed 
that a moderately massive disk (10~ 5 — 10~ 2 M Q ) can 
cause the precession of the central compact star. Trumper 
et al. (1983) used the free precession of the neutron star 
to explain the 35— day period. Qiao & Cheng (1989) sug- 
gested that instead of the free precession of the neutron 
star, the periodic variation is due to the precession of the 
accretion disk driven by a normal star. 
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In this paper, we present an accretion disk model to 
account for the precession of pulsars. The total torque ex- 
erted on the pulsar by the fossil disk, which causes the 
pulsar to precess, is nonzero, due to the asymmetry in 
pulsar shape and nonalignmcnt between the disk rotation 
axis and the axis of pulsar moment of inertia. Hence, we 
present this forced precession model (disk— induced pre- 
cession) , other than the free precession model. 

In fact, such a system, i.e. a pulsar with an accretion 
disk, has been invoked to describe the nature of two kinds 
of isolated X— ray pulsars: the anomalous X— ray pulsars 
(AXPs, which have persistent X— ray pulsations but few 
burst) and the soft gamma— ray repeaters (SGRs, which 
have multiple bursts of soft gamma— rays) (see Xue et al. 
2003 for a review). The most remarkable difference be- 
tween these two kinds of radio— quiet X— ray pulsars and 
the radio pulsars is that their X— ray luminosities (usually 
Lx ~ 10 35 ergs/s) are several orders of magnitude higher 
than their rotational energy loss rates. To explain this and 
other exotic facts, two alternative models, viz. the mag- 
netar models (in which the pulsar is suggested to be a 
neutron star with super— strong magnetic field) and the 
accretion disk models, were proposed (see Xue et al. 2003 
for a review). Under the framework of the accretion disk 
model, a problem would come if the precession of radio 
pulsars (e.g. PSR B1811-11 and B1642-03) was really 
caused by the accretion disk: why do some pulsars with 
disks behave as ordinary radio pulsars while others behave 
as AXPs/SGRs? In this paper, we discuss this problem in 
section 3 after presenting the detailed calculation in sec- 
tion 2. 



2 



Please give a shorter version with: \authorrunning and/or \titilerunning prior to \maketitle 



2. Pulsar precession caused by an accretion disk 



For simplification, we assume c ~ d « R, so that 



2.1. Geometry and parameters 

Consider the neutron star as a rotation ellipsoid with the 
principal moment of inertia ii = I 2 ^ Is and the corre- 
sponding radii a = a =/= b. The geometry of the system 
is shown in Fig.l: O—xyz is the coordinate system fixed 
with the neutron star, Ox, Oy and Oz arc the principal 
axes of moment of inertia of the neutron star, the axis 
OA is the rotation axis of the neutron star (for simplic- 
ity, we assume the axis Oz is aligned with the axis OA in 
this paper; in this case, free precession doesn't exist and 
we can see how the disk makes the neutron star precess 
more clearly), 0—^r]( is the coordinate system fixed with 
the disk and the disk plane is O^rj, 9, tp and cj> are Euler 
angles, uj is the angular velocity of neutron star rotation. 
We will derive the formulae of the gravitational potential 
in the neutron star— disk system and the angular veloc- 
ity of neutron star precession. And then, as an example, 
we model the precession period of PSR B1818— 11. The 
parameters about PSR B1828— 11 are as follows: the ro- 
tational period of this pulsar P = 0.405s (w = 2-jt/P), the 
period derivative P — 6.0 x 10~ 14 ss _1 , the precession pe- 
riod P prc w 1000 days, the magnetic field B = 5.0 x 10 12 G 
(Stairs et al. 2000). 



2.2. Equations and numerical results 

The gravitational potential between the neutron star and 
the disk is 
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where M, a and b are the mass, the equatorial and longi- 
tudinal radii of the neutron star, respectively, M is the 
mass of the disk, c and d are the inner and outer radii of 
the disk, respectively. 

The moment of inertia of the neutron star is I z — 
(8/15)7rpa 4 6, where p — iM/ (Aira 2 b) is the density of the 
neutron star. Since the rotation angular velocity of the 
neutron star ui is much greater than the precession angu- 
lar velocity of the neutron star </>, the angular momentum 
of the neutron star can be written as Lq ~ I z lu. According 
to the Lagrange's equation L <fism9 = —dV/d0 (Qiao & 
Cheng 1989), we get the precession angular velocity of the 
neutron star: 



3G M cos 9 b 2 
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In terms of the study on the stability of the Maclaurin 
spheroids, we can obtain the oblateness of this pulsar: 
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c = 



= 2.27 x 10" 
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Fig. 2 shows the relations between M and R derived 
from the 1000— day precession period for a group of 9. 
In principle, Mq may vary from less than O.OOIMq to 
even 0.1M Q (The amount of the fall— back material (Mfb) 
which forms the disk after supernova explosion is uncer- 
tain, however, Mfb < 0.1M© is supposed to be reasonable, 
see Chatterjee et al. 2000), and R from about li? co (or 
smaller) to 2 ~ 3i? co when varies within the range of 
~ 90°, where R co is the corotation radius of the disk 
with the pulsar. Thus, the precession period of the PSR 
B1828-11 is well elucidated. 



3. Possible evolutionary scenario in the 
pulsar-disk system 

One may ask where the accretion disk comes from. 
Usually, a disk around the pulsar can be formed: 1. by 
accreting the matter from a companion star in the binary 
system (approach I); 2. by the fall— back of the material 
after a supernova explosion (approach II). In the former 
process, there are two ways of mass transfer: matter flow- 
ing through the first Lagrange point due to the critical 
Roche Lobe abounding with the companion matter, or by 
the stellar wind (Nagase 1989). 

An interesting issue is that what the pulsar in the 
pulsar— disk system will behave as, radio pulsar with pre- 
cession or AXP/SGR? In the accretion model of X— ray 
binary, accreting process may experience some "phases", 
e.g. "radio pulsar phase", "sleeping phase" , "propeller 
phase" and "accretion phase" (Illarionov & Sunyaev 1975, 
Henriches 1983). The case is suggested to be similar for a 
pulsar— disk system formed through either approach I or 
II. The PSR B1828-11, B1642-03, AXPs and SGRs are 
all isolated "neutron stars". What are the links between 
them? Here we present a possible series of evolutionary 
tracks for the pulsar— disk system (cf. Chatterjee et al. 
2000 and Henriches 1983). 

(a) After the supernova explosion, the newborn neu- 
tron star rotates so fast that its magnetospheric radius 
R m is much larger than its light cylinder radius Rlc- The 
coupling between the disk and the neutron star is unlikely. 
The young pulsar initially emits a very strong relativistic 
wind together with highly energetic magnetic dipole radi- 
ation. The pulsar wind and photons blow a cavity in the 
fall— back flow and no matter can accrete, as is similar 
in the binary system (Henriches 1983). Thus the neutron 
star may live in the "radio pulsar phase", i.e. the radio 
emission can propagate away. 

(b) The neutron star slows down as its rotational en- 
ergy reduces. Its radiation and the high energy particle 
flow become weakened, thus the cavity decreases and the 
radio emission is finally quenched. When the "radio pulsar 
phase" ends, accretion down to the surface of the neutron 
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star is not expected immediately. The so— called "sleep- 
ing phase" takes place. This phase lasts until the pulsar 
rotation slows down enough to make accretion possible. 

(c) Subsequently, the neutron star enters the "pro- 
peller phase" since accretion onto the neutron star is 
inefficient due to centrifugal forces acting on the mat- 
ter. As the rotation slows down further, accretion will 
be permitted and the neutron star will behave as an 
accretion— driven X— ray source. Gradually, the system en- 
ters a quasi— equilibrium "tracking phase" (happening at 
the time of ttrans), in which the spin of the star roughly 
matches the rotation of the disk at R m . In this phase, 
the neutron star will be X— ray bright due to accretion. It 
is possible that the advection— dominated accretion flow 
(ADAF) will occur at the time of t adaf i when much of 
the mass will be ejected prior to reaching the neutron star 
surface. The X— ray luminosity of the neutron star will de- 
cline rapidly with time due to Mx (the rate of accretion 
onto the star surface) diminishing. Therefore, the bright 
AXPs may exist only between ttrans and 2£adaf (about 
5000 - 40000 years). 

(d) Finally, accretion ceases since the matter of the 
disk is exhausted. Thus no accretion— driven X— ray pulsar 
can be observed and neither can radio emission be emit- 
ted owing to the long rotational period of the neutron star 
(Qiao & Zhang 1996, Zhang et al. 2000). Consequently, 
the neutron star may be observed as a dim isolated neu- 
tron star, a weak soft X— ray source powered by rotational 
energy loss or thermal radiation. 

4. Conclusions and discussions 

The accretion disk model for periodic timing variations of 
pulsars (e.g. PSR B1828— 11) can elucidate the observed 
precession period. According to the analysis of the possible 
evolutionary scenario discussed above, the PSR B1828— 11 
may be just in the "radio pulsar phase" before the "sleep- 
ing phase" and the "propeller phase" . If our model proves 
to be true, the link between radio pulsars and AXPs may 
be established and the fossil disk model for AXPs may be 
strengthened. It also suggests that the neutron star— fossil 
disk system might exist. Some discussions are listed below. 

1. Sarazin et al. (1980) pointed out when the angu- 
lar momentum of the disk is larger than that of the neu- 
tron star, the neutron star will precess due to the effect 
of the accretion disk and the disk will maintain its orien- 
tation. This is just our case for the precession of the PSR 
B1828-11. 

2. Why is the radio emission expected in our model 
but not expected in the case of close X— ray binary sys- 
tem? Illarionov & Sunyaev (1975) suggested that the large 
optical depth of the stellar wind for free— free absorption 
hinders the detection of radio pulsars. In our model, a 
fossil disk is formed and there is no enough stellar wind 
to absorb the radio emission. A case in point is the PSR 
B1259— 63, which has a B c companion star with high el- 
lipticity and long orbital period. This neutron star is ob- 
served as a radio pulsar when the companion encircles far 



away from it. But no radio emission is detected when the 
companion approaches very near the pulsar (Johnston et 
al. 1992). 

3. The condition when the pulsar radio emission 
quenches is not clear yet. A possible case may be that 
the cavity radius becomes smaller than the light cylinder 
radius (Henriches 1983). In our model, although the coro- 
tation radius of the PSR B1828-11 R co = (GM/w 2 ) 1/3 w 
9.19 x 10 7 cm is one order of magnitude smaller than its 
light cylinder radius i?LC = c/u = 1.93 x 10 9 cm, and the 
disk radius R approaches 2 — 3i? co , this neutron star is still 
in the "radio pulsar phase" . This presents an example for 
further investigations. 

4. We have ignored the angular momentum of the 
pinned supcrfluid in the dynamic equations of the calcula- 
tion. Can the pinning force be strong enough to damp any 
precession of neutron stars? In principle, detailed compu- 
tations may determine whether the vortices are pinned, 
but a definitive conclusion on the nature has not been 
reached yet due to various uncertainties in the microscopic 
physics. Based on the calculation of Shaham (1977) and 
Sedrakian et al. (1999), Stairs et al. (2000) conclude that 
the torque— free precession of PSR B1828— 11 should be 
damped out due to the dissipation caused by the vortices 
pinning to the stellar crust on the timescale of several hun- 
dred precession periods. Nevertheless, a further consider- 
ation of Shaham (1986) indicates an unpinned supcrfluid 
in the case of precession. It is also noted by Link & Cutler 
(2002), that hydrodynamic forces present in a precessing 
star are probably sufficient to unpin all of the vortices of 
the inner crust. In addition, the 35— d periodicity of the 
accreting binary system (Her X— 1) has been interpreted 
as precession (Shakura, Postnov & Prokhorov 1998, and 
the references therein), which could be another example of 
forced precessions. We have therefore a tendency to sug- 
gest that, the observation of PSR B1828— 11 precession 
may imply that the vertex pinning in this pulsar is much 
weaker than that predicted previously, at least it is weaker 
than the possible disk— driven precession effect. 

5. In fact, the phenomenon of the secular periodic tim- 
ing variations, which also exists in many other radio pul- 
sars (e.g. Crab and Vela pulsar), has been investigated 
both observationally and theoretically (Sedrakian et al. 
1999 and the references therein). If a neutron star were 
a rigid or semirigid solid star, it would precess (Pines & 
Shaham 1974); as is mentioned above, the interaction be- 
tween the superfluid vortices and the stellar crust needs 
further considerations. Ruderman (2001) pointed out that 
this is a crucial problem which will raise doubts about 
canonical descriptions of neutron star structure until it is 
solved. No matter the precession is torque— free or forced, 
the neutron star structure still needs to be investigated 
more. Instead of neutron star, Xu (2003) proposed a solid 
strange star model, which may be possible and also needs 
to be investigated more. 

6. Our model cannot exclude the possibility of a com- 
panion star (instead of a fossil disk) which can produce 
the precession phenomenon as well. In this case, two limits 
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should be considered: the Roche Limit d,R (if the compan- 
ion approached the neutron star closer than this limit, it 
would be torn asunder by tidal forces) and the Instability 
Limit d\ (if the company strayed farther than this limit, 
it would escape due to the differential perturbations from 
other celestial bodies). The details of this case will be stud- 
ied in another paper. 
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Fig. 1. The geometry and parameters for the calculations. 
In our calculation, the rotation axis OA is aligned with the 
principal axis Oz of moment of inertia of the neutron star. 



0.5 




1 1.2 1.4 1.6 1.8 2 



R/R 

CO 

Fig. 2. The relations between Mo and R derived from the 
1000— day precession period for a group of 9. 



